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Abstract 


We present the quantum equation and synthesize an optimal control proce- 
dure for this equation. We develop a theoretical method for the analysis of 
quantum optimal control system given by the time depending Schrédinger 
equation. The Legendre wavelet method is proposed for solving this prob- 
lem. This can be used as an efficient and accurate computational method 
for obtaining numerical solutions of different quantum optimal control 
problems. The distinguishing feature of this paper is that it makes the 
method, previously used to solve non-quantum control equations based 
on Legendre wavelets, usable by using a change of variables for quantum 
control equations. 
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1 Introduction 


With the advent of the twentieth century, the inability of classical physics in 
the fields of relativity and microscopy led to the outburst of quantum physics. 
After the development of quantum physics, the issue of quantum control 
was inspired by experimental advances and issues raised in sciences such as 


*Corresponding author 

Received 19 October 2020; revised 5 May 2021; accepted 30 May 2021 

Mojtaba Rahimi 

Department of Mathematical Sciences, Yazd University, Yazd, Iran. e-mail: rahimimo- 
jtaba1358Q@gmail.com 


Seyed mehdi Karbassi 
Department of Mathematical Sciences, Yazd University, Yazd, Iran. e-mail: smkar- 
bassiQ@yazd.ac.ir 


Mohammad Reza Hooshmandasl 
Department of Computer Science, University of Mohaghegh Ardabili, Ardabil, Iran. 
hooshmandasl@uma.ac.ir 


333 


334 Rahimi, Karbassi and Hooshmandasl 


quantum chemistry, quantum optics, quantum information, and atomic and 
molecular physics [7, 10, 24, 27, 34]. Recently, classical control methods such 
as optimal control, robust control, Lyapunov control, and feedback control 
for quantum systems have been studied and expanded [5, 16, 19, 20, 22, 25]. 
In quantum control, the main goal is to effectively control the system from 
an initial state to a desired final state using external control fields. However, 
control in quantum systems is at the beginning of the road, and further 
research is needed. In particular, optimal control in quantum systems is of 
particular importance as one of the most widely used issues [3, 9, 12, 15, 17, 
23, 26, 31]. 

Along with the development of analytical control methods in quantum 
systems, numerical methods have also been developed. Numerical methods 
that have been considered in solving classical equations and problems are also 
studied [8, 11, 18, 28]. One of the most useful numerical methods for solving 
differential problems is wavelets-based numerical method [2, 14, 21, 29, 33]. 
This article tries to provide a useful way to solve optimal quantum control 
based on wavelets method. 

Let H be a finite- or infinite-dimensional Hilbert space of a quantum system 
and let W denote the state of this system. Then the Schrédinger equation 
can be found as follows [13, 32]: 


(1) 


where W € H is the state variable, h is the Planck constant, and H isa 
self-adjoint Hamiltonian operator in H. 

In every physical system, energy is an important quantity. In quan- 
tum systems, Hamiltonian H is corresponding to energy, then we can write 
He, = Eye;,i = 1,...,N, where F; are eigenvalues, e; are eigenvectors of 
the physical system concerned, and N is the dimension of H. In this paper, 
we suppose that e; is an orthogonal basis; then we expand a state vector V 


as follows: 
N 
i=l 


The population of energy states of level i are the quantities |q); 
When a quantum system is operated by an external field, the Schrédinger 
equation in (1) is modified as 


|?. 


OW 
ih = (Ho + e Hyun (t))¥, 


where Hy : H+ H, which is the internal Hamiltonian, and the Hamiltonian 
linear operator Hy, :H+> H describes the coupling of the system to external 
fields u,(t). In this paper, for simplicity, we consider a quantum system with 
one control, and set k = 1. Then we can write 
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Free! = (Ho + Ayu(t))Wv. 
Ot 
The design of controls in quantum systems is considered for energy- 
efficient population transfer. These controls perform the desired transmission 
and in addition optimize a specific performance index. The specific energy 
performance index considered in this section is shown below: 


4 
min J[u] -f u?(t)dt. 


This cost function for control has been widely used in the literature on optimal 
control of quantum systems as a part of various objective functions. It is a 
measure of the energy expended to create a control field. 


The main object of this paper is to present an efficient numerical algorithm 
based on the Legendre wavelets methods to solve the following optimal control 
problems of the form: 


T 
min J[u] a u?(t)dt, 
0 
subject to the dynamical quantum systems 


iW = (Hy +Vult))W. 


The initial condition for the above equation is 


It is worthy to note that V = Aj. 


2 The hat function 


In this section, we introduce a family of basic functions, namely, the hat 
functions. An n-set of the hat functions is defined on the interval [0,7] as 
follows: 


GAVE (fj -l)k <t< jk, 
hj (t) as GtD ket jk) < t < (J + 1)k, 
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pS). Tuapege rT 
init) ={, ° O.wW oe 


where k = —'. For hat functions, we can write hj(jk) = 6;;, where 6 is the 
Kronecker delta. By the definition of the hat functions, we can expand any 


function like g(t) € L?[0,T] as follows: 
n-1 
g(t) = D0 gyhy(t) = GT H(t) = H7()G, (2) 
j=0 


where 
G = [90, 91; aes »In—-1] 
and 
H(t) © [ho(t), hi(t),..., hn—1(t)]*- 


When we use the hat functions for the g(t), it can be observed that 


Now, we introduce another family of basic functions, namely, Legendre 
wavelet functions. The set of these functions is an orthogonal set on the 
interval [0,1] with the weight function w(t) = 1. If P,(t) is the Legendre 
polynomials of degree k that are orthogonal on the interval [—1, 1] with re- 
spect to the weight function w(t) = 1, then we can write Legendre wavelets 
as follows [2, 14]: 


0, O.W. 


V2k + 12? P,(2"*'t- 214-1), t€ (EA, Al, 
w(t) = 


In fact, wiz(t) = w(l,k,n,t), where | = 1,2,...,2" and n is an arbitrary 
positive integer. 

For any arbitrary function, like g(t) defined over [0,1] and square-integrable 
over [0, 1], we can expand g(t) by the Legendre wavelets as follows: 


Co CoO 
g(t) _ S- > ae Wig (t). 
I=1 k=0 
By approximating the above infinite series, we can write 


2” K-1 
g(t) ~ Qik Wik (t) = A Wit), (5) 


where A and W(t) are k = 2"K column vectors. 
For the index Jk, we can write 7 = k(I—1)+k+1. Then ay, = a; and 
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Wik = w;. Thus (5) can be written as 


k 
g(t) = S) ajyw;(t) = ATW(t), 
l=1 


where 
AS [a1,a2,...,a,]” 
and 
W(t) = [wi(t), wo(t),-.-, wet). (6) 
Now by taking t; = naa as the collocation points into (6), we can write 
re 1 
Peg = I) ,W(1)] 


It can be simplify verified that the Legendre wavelets can be expanded in k 
in terms of the hat function by using (2) and (3) as follows: 


W(t) ~ Py gH (0) (7) 
Theorem 1. Suppose that f(t) ~ F7 H(t) and that g(t) ~ G7 H(t). Then 

f(t)g(t) ~ STH), (8) 
where 5;; = (F.G);; = FijGi; denotes pointwise product of F and G. 


Proof. By applying (2) and (3) for f(t) and g(t), we can write 


n—-1 k-1 
g(t) ~ GPH(t) = > gjhy(t) = © gGk)hy(t), 
j=0 j=0 


n-1 k-1 


f(t) = FT H(t) = D0 fyhy(t) = Do FGR)HG(t)- 
j=0 


7= 


Then by using the point wise product, we have 


k-1 


F()g(t) = D0 fk) g(Gk)hy(t) = DT H(t), 


j=0 


which completes the proof. 


Corollary 1. Suppose that g(t) ~ G7 H(t) by hat functions. Then, for any 
integer number n > 2, we have 


(9(t))” = [90 91.---.9¢_,)H@). (9) 
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Proof. For n = 2 by Theorem 1, we have 

(a(t)? © (90; 7>-=+382_,]2@; 
then by induction for n > 2, we have 

(9())” = [90 91,---.9¢_,JH@), 


which completes the proof. 


Theorem 2. Suppose that g(t) ~ A? W(t) and f(t) ~ B?W(t) by Legendre 
wavelets. Then we can write 


f(t)g(t) ~ QT PW (8), (10) 


where Af = AP, ;,, BT = BTP; 


kx ko and Q = A,.By. 


xk? 
Proof. By Theorem 1 and equation (7), we have 


g(t) ~ ATW(t) ~ ATP, ; H(t) = Af H(d), 


f(t) = BY W(t) ~ BYP, H(t) = Br H(t), 
and then 


f(t)g(t) = (Ar-Bi)? H(t) = QTH(t) = QTPLLW(t). 


In Theorem 1, the multiplication of two functions is obtained according to 
the hat functions, and in Theorem 2 by using Theorem 1, the multiplication 
of two functions is obtained according to the Legendre wavelets. If g(t) ~ 
A™W/(t) by the Legendre wavelets, then by Theorem 2 and Corollary 1, we 
can write 

(g(t))” = [af 43, ...,a2]P-", W(t). (11) 


3 Analysis of the proposed method 


In previous works, using Legendre wavelet, numerical solutions for non- 
quantum control equations have been obtained. In this section, we try to 
implement one of these methods, which is based on the Legendre wavelet, 
for quantum control equations. To overcome this problem, first change of 
variables to make the equation usable is performed. 

In this section, we consider the quantum control systems of the form 
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e 
min J[u] =| u?(t)dt, 


subject to the dynamical system 


i = (Hp +Vult))¥, (12) 


with the initial condition U(0) = Wo. 
First, we introduce a change of variables that generalizes 


v= e Hoty. 
such that for (12), we obtain 
4(—iHye “2a + e-otg) = (Hy + Vult))e "es Hype te + oot 
= Hye Pete + Vu(t)eFo'g. 


Thus we can write 
& = eb Hotyy (te Hoty. 


By considering E(t) = e’40*Ve-*#0!, equation (12) can be written as 
& = E(t)u(t)ax(t), (13) 


and the new initial condition is 7(0) = Wo. 
By the Legendre wavelets for the derivative of the state variable ¢ and the 
control variable u(t), we can write 


&~ XTW(t) (14) 
and 
u(t) ~ UT W(t), (15) 
where 
Ur = [u1, U2, , Up] 
and 
XT = [e1,42,... ytzl. 


We can apply integration on both sides of (14) and considering the initial 
condition 
a(t) ~ XTQW(t) + Wo. 


Let 7 be the coefficients vector of the unit function. Then 
x(t) = (X7Q + Von" )W(t). (16) 


Now (16) can be written as 
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a(t) ~ ATW(t) = [a1,a2,...,a;]W(t). (17) 


In this part, by applying (11) and (9) for the above approximation, we have 


(w(t))” = [G7 4g,...,@2]P-", W(t) = [a7 ag,..., a2]H(), 
and also for (15), we have 
(u(t))” = (a, a... 02] PoL W(t) ~ [at a... aH). (18) 
We can also approximate E(t) by Legendre wavelets as 
E(t) ~ E’ W(t), (19) 


where 
ET = [e1,e2,.. ‘ye |: 


In this part, by using (8) and (10) for E(t)u(t)x(t), we have 
E(t)u(t)a(t) ~ [6c t1, éotinda,... , Egtig%,]P-", W(t) 
 [eyuid1, €otioe,++-,eptig Be) A (4) 
= ATH), (20) 


also by using (18) for n = 2, we have 


(u(t)? ~ [&7, @,...,@2)Po'. W(t) & (a7, 3,..., JH) = AZ A(E). (21) 


Now by using (21), the index J can be written as 


J ~ J[U] = AFQ, (22) 


= if rotbat, f° hi(t —e hg_,(t)dé]. (23) 


By applying (20) for (13) and (12 ), we can write 


where 


XT — A, Pol. ~0. (24) 


In this part, by Lagrange multiplier method for minimization index J in (22) 
subject to systems of algebraic equation (24), we can write 


J[X,U, E] = J[U] + X7 - APL =AZQ+XT-AyPo.L, (25) 


where 
Ie Diy dase cig he) 
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Hence L is the vector of Lagrange multiplier. 
For minimizing by the Lagrange method, the necessary conditions are 


OJ 

ax ~% 

OJ 

OF 26 
ary = (26) 
OJ 

ao 


By the Newton iteration method, we can solve equations of (26) for X, U, 
and L. Then the approximation of z(t) and u(t) can be determined by (17) 
and (15). 


4 Proposed algorithm 


The object of this algorithm is designed to solve the Schrédinger equation: 
Input: T (final time), N( dimension of H), Ho, and Hy . 

Step1: Make a change of variable V = e~*/0tx in dynamical system (12) to 
obtain (13). 

Step2: Define x(t), u(t), and E(t) by (17),(15), and (19), respectively. 
Step3: Write the index J as (22) and (23). 

Step4: Compute dynamical system (13) by applying (20) to the form (24). 
Step5: Compose new index J[X,U, L] as (25) by (24) and (23). 

Step6: Solve equation systems in (26) and obtain X and U. 

Step7: Compute x(t) and u(t) by (17) and (15). 

Step8: Compute VU = e~*oty, 

Output: The approximate solution u(t) and W. 


5 Numerical experiments 


Example 1. In this example, we consider the two-level system iv = (Ho + 
Vu(t))W, where W € C? as follows: 


_ EF, 0 _ 0 vie 
Ho (S im) Ya (4,0) 


The concept of optimal control of two-level quantum systems was presented 
in [1, 3, 4, 6, 8, 11, 18, 30]. In the most of these researches, optimal control 
is constructed on the basis of geometric arguments. If we suppose EF, = 2, 
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Eg = —2, and v12 = 2+3:2, then we can obtain the numerical results as follows: 


Table 1: Approximate and exact value of u(t) in Example 1 for different values of t. 


t approximate of u(t) exact of u(t) error 
1 -0.0359 -0.0459 0.0120 
5 -0.0496 -0.0578 0.0102 
10 0.1010 0.0999 0.0051 
15 -0.0922 -0.0878 0.0016 
20 0.0663 0.0594 0.0119 
25 0.0286 0.0066 0.0220 
30 -0.0630 -0.0635 0.0020 


Table 2: Approximate and exact value of a(t) in Example 1 for different values of t. 


t approximate of y(t) exact of y(t) approximate of we(t) exact of we(t) error of q(t) error of y(t) 


ail 1.0097 1.0078 0.0690 0.0590 0.0109 0.0190 
5 0.9869 0.9769 0.2664 0.2564 0.0108 0.0190 
10 0.8666 0.8766 0.4681 0.4764 0.0110 0.0113 
15 0.7197 0.7297 0.6656 0.6756 0.0120 0.0160 
20 0.5513 0.5413 0.8645 0.8445 0.0110 0.0230 
25 0.3033 0.3133 0.9270 0.9470 0.0120 0.0220 
30 0.0830 0.0737 1.0105 1.0005 0.0123 0.0130 
0.15 T 


Exact and approximate solutions 


0 5 10 15 20 25 30 35 


Figure 1: Plots of approximate and exact results of control variable for Example 1 
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Exact and approximate solutions for x(t) 


35 


Figure 2: Plots of the approximate and exact results of state variable for Example 1 


io 


= 
fo) 


= 
ns 


-_ 
ip 


Oo 


Exact and approximate solutions of population 
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Figure 3: Plots of the approximate and exact results of population for Example 1 


We solved the above problem by our proposed algorithm with n = 2 and 
K = 10 or k = 40. The stopping condition is |u*(t) — u(t)| < 5 x 1072, 
where u*(¢) and u(t) are approximate and exact results of control variable, 
respectively, and by the stopping condition, we have 73 iterations. The plots 
of the approximate and exact values of u(t) is shown in Figure 1. The plots 
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of the approximate and exact values of ~(¢) and 72(t) are shown in Figure 
2. The plots of the approximate and exact values of |7,(t)|? and |2(t)|? are 
shown in Figure 3. 

Approximation and exact results of control variable u(t) and their absolute 
errors for different t are presented in Table 1. In fact, because the index J(t) 
in each step depends on the control variable, so by approximating the control 
variable in each step, the index J(t) is also approximated. Approximation 
and exact results of state variables ~,(t) and w2(t) and their absolute errors 
for different t are presented in Table 2. The above results show that using 
the algorithm and method mentioned numerically is very useful and efficient. 
The main advantage of this method is to provide a simple solution based on 
classical numerical methods in the field of optimal quantum control. In [32], 
an index similar to the index used in this article has been used and numerical 
and graphical results have been obtained. Carefully in these results, it is 
observed that the solutions obtained in this article have been obtained with 
less repetition and more accuracy. 


Example 2. In this example, we consider the three-level system iV = (Ho + 
Vu(t))W, where UW € C? as follows [32]: 


FE, 0 0 0 vy2 013 
Ho = 0 Eg, 0 ; V= Uio 0 U23 
0 0 Es Viz 93 0 


The concept of optimal control of three-level quantum systems were presented 
in [3, 4, 13, 32]. In the most of these works, the optimal control is constructed 
on the basis of geometric arguments. If we suppose EF, = 2, Ey = 0, F3 = 6, 
Vyg =t4+1, vig = 4, and vo3 = 2+ 3i, then we can obtain the numerical 
results as follows: 


Table 3: Approximate and exact values of u(t) in Example 2 for different values of t. 


t approximate of u(t) exact values of u(t) error 

2 0.0074 0.0074 3.3 x 10-4 
10 0.0029 0.0029 5.2 x 1074 
20 -0.0107 -0.0105 2.0 x 10-4 
30 -0.0002 -0.0001 1.1 x 1074 
40 0.0062 0.0062 3.3 x 1074 
50 -0.0032 -0.0032 2.3 x 10-4 


60 0.0060 0.0063 4.1 x 10-4 
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Table 4: Approximate and exact values of U in Example 2 for different values of t. 


t approximate of approximate of approximate of error of error of — error of 
wr(t) wa(t) wa(t) vr (t) w2(t) w(t) 
2 1.0188 -0.1011 0.0731 0.0200 0.0220 0.0220 
10 0.9837 -0.3944 0.2748 0.0140 0.0154 0.0220 
20 0.9016 -0.6478 0.5012 0.0210 0.0141 0.0123 
30 0.7168 -0.7977 0.6782 0.0220 0.0215 0.0144 
40 0.5713 -0.6806 0.8288 0.0180 0.0125 0.0214 
50 0.3356 -0.4449 0.9638 0.0000 0.0112 0.0124 
60 0.0799 -0.0525 0.9767 0.0225 0.0114 0.0126 
0.015 1 


0.01 


0.005 


-0.005 


Exact and approximate solutions of u(t) 
lo} 


-0.01 


-0.015 
0 


70 


Figure 4: Plots of approximate and exact results of control variable for Example 2 
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Exact and approximate solutions for state variable 


70 


Figure 5: Plots of the approximate and exact results of state variable for Example 2 
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Figure 6: Plots of the approximate and exact results of population for Example 2 


We solved the above problem by our proposed algorithm with n = 2 and 
K = 10 or k = 40. The stopping condition is |u*(t) — u(t)| < 5 x 1074, 
where u*(t) and u(t) are approximate and exact results of control variable, 
respectively, and by the stopping condition, we have 97 iterations. The plot 
of the approximate and exact values of u(t) is shown in Figure 4. The plots 


Solving quantum optimal control problems by wavelets method 347 


of the approximate and exact values of y(t), We(t), and w3(t) are shown in 
Figure 5. The plots of the approximate and exact values of |q1(t)|?, |¢2(t)|? 
and |w3(t)|? are shown in Figure 6. 

Approximation and exact results of control variable u(t) and their absolute 
errors for different t are presented in Table 3. Approximation and exact 
results of state variables ~1(t), we(t), and w3(t) and their absolute errors 
for different t are presented in Table 4. Figures and tables obtained above 
show that the numerical method used in this paper is simpler and more 
useful than other methods. In [32], an index similar to the index used in 
this article has been used, and numerical and graphical results have been 
obtained. Investigating these results carefully, it can be observed that the 
solutions obtained in this article endure less repetition while having more 
accuracy. 


6 Conclusion 


Today, the issue of quantum optimal control is one of the most widely used 
issues in many basic sciences and engineering. At the same time, numerical 
methods are one of the most useful solutions to these problems. In this paper, 
a numerical method based on wavelets was proposed to solve the problem of 
optimal quantum control. This method was benefited by using topics related 
to applied mathematics in the field of classical equations and presented usable 
in the field of quantum systems. The above results showed that using the 
algorithm and method mentioned numerically is very useful and efficient. 
The merit of this method is to provide a simple solution based on classical 
numerical methods in the field of optimal quantum control. 
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